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ABSTRACT 

The phenomenon of clustering of galaxies on the basis of correlation functions in an ex- 
panding universe is studied by using equation of state, taking gravitational interaction 
between galaxies of extended nature into consideration. The partial differential equa- 
tion for t he extended mass structures of a two-point correlation function developed 
earlier by llabal et al.l (|2006[ ) is studied on the basis of assigned boundary conditions. 
The solution for the correlation function for extended structures satisfies the basic 
boundary conditions, which seem to be sufficient for understanding the phenomena 
and provides a new insight into the gravitational clustering problem for extended mass 
structures. 

Key words: gravitational clustering, thermodynamics, extended structures, correla- 
tion function, structure of universe. 
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1 INTRODUCTION 

The universe is dominated by matter in gravitational in- 
teraction into the spatial configuration consisting of galax- 
ies, g roup of galaxies, super-clu sters and even larg er struc- 
tures [Shane&iWkiinei (119671 ) and lMaddo:3 (Il990l '). One of 
the standard ways in which structures in the universe were 
formed (i.e; the departure from randomness and homogene- 
ity) is described by means of correlation functions. From 
observation the galaxy-to-galaxy correlation function (^gai is 
known to scale as ^gai = r~'' , where the exponent 7 ranges 
from 1.6 to 1.8. The calculation of 7 tacitly assumes the 
universe to have evolved to the stage, where the initial pri- 
mordial matter has been formed into the observed galaxies 
and these galaxies are coupled to the expansion of the uni- 
verse. We take our clue by analogy with the well established 
theory of interacting gases iGoodsteinI (|l985l ). where the rel- 
ative arrangements of the atoms and molecules (galaxies in 
our case) are accurately described in terms of the princi- 
ples and methods of thermodynamics. We apply the same 
techniques to a system made up of many galaxies in grav- 
itational interaction, subject to fluctuations emerging from 
the intrinsic properties of the gravitational interaction. A 
number of theories of the cosmological many body prob- 
lem have be en developed mainly from t he thermodyn amic 
point of viewlSaslaw fc HamihonI ligS^ I. ISaslawl (|200(]| 'l and 
llgbal et al.l (|2006l ) . The theory makes use of the equations of 
state along with the correlation functions for the extended 
mass structures for the development of a semi-analytical 
model. This can even be done by solving a system of Li- 
ouville's equation or BBGKY-hierarchy equations and have 
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been discussed by iPeeblej (|l980l ) and ISaslaw fc HamiltonI 
("1984). But BBGKY-hierarchy equations are too compli- 
cated to handle for higher order correlation functions. We 
extract the possible information about lowest order i.e; two- 
point correlation function for galaxies with real extended 
mass structures clustering gravitationally in an expanding 
universe. Although, th e same information about ^2 have ear- 
lier been developed bv llabal et al. but here we have 
studied and gather the related information about ^2 for ex- 
tended mass structures. It is important to note that galax- 
ies with point mass consideration is only an approximation. 
In fact, galaxies have real extended structures, where dark 
matter especially is having an important contribution. 



2 TWO-POINT CORRELATION FUNCTION 
FOR EXTENDED MASS STRUCTURES 

The general pair of equations of state for internal energy 
(Ue) and pressure (Pe) for extended mass structures is given 
as; 



Ue = -NT 



Pe = 



NT 



2TvGm^N^ 
V 

2TTGm'^N'^ 



e(n,r,r) 1 + 



dV 
Attv 



(2) 



V 31/2 

Equat ions( l) and (2) r e presen t the equations of state iHilll 
(|l956l ) and llabal et al.l (|2006l ). The measuring correlation 
parameter (be) for extended mass structures is given by 



2nGm^N 
3VT 



(3) 
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Here n = ^ is the average number density of extended 
mass particles (galaxies) with t as softening parameter and 
is taken between 0.01 to 0.05 (in the units of total radius), 
T is the temperature, V the volume, G is the gravitational 
constant, ^(n, T, r) is the two-point correlation function for 
extended mass structures and 'r' the intergalactic distance. 
It should be noted that these expressions assume a large 
volume V for their validity. 

The Maxwell's thermodynamic equation relating inter- 
nal energy (Ue) and pressure (Pe) for extended mass struc- 
tures is given by 

(9Jk\ ^t(^] -P 

\dVjT,N \dTJv,N 

From equations (1) and (2), we have 
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with Ci as constant. 

Similarly, in order to find the solution of the equation (11) 
we can proceed as; 



+ f r dR 
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Equation (13) can be correct only, if both sides of it are 
equal to the same constant say ' Z'j^, 

+ r-^ ( r dR 



\R dr 

TdT 

— Zjf\J 
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The possible solutions of equations (14) and (15) are; 

R{r) = C2r^ {e^ + r^)^ 
r(r) = C:iT^" 



(14) 
(15) 

(16) 
(17) 



with C2 and C3 as constants. 

With the help of the above possible solutions, the two point 
correlation function defined by equation (8) is 



6 {n, T, r) = C{n) 3 T^" (e' + r") — 
where C = C1C2C3 



(18) 
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Using the equations (5) and (6) in equation (4), and 
then differentiating the required equation with respect to V, 
the two-point correlation differe ntial equation for e xtended 
mass structure takes the form as llabal et al.l (|2006l ): 



Sn—- + T 
an 



J.2 _|_ 



9e 



96 „ 

r — 

dT dr 



(7) 



Here ^2 depends on the variables n, T, and r. In order 
to look for the possible solution of this equation, we write 
^2 as product of three variables defined by; 
e(n), r(r), and R(r). 



^2{n,T,r) = e{n)r{T)R{r) 



(8) 



Calculating the required terms from equation (8), equa- 
tion (7) gets modified to 



3n de(n) r dR{r) 



dr(T) 



Q dn R dr T e'^ + r'^ dT 

As both sides of equation (9) are functions of different 
variables, so equation (9) will be true only, if both sides can 
be equal to the same constant say 'Z'. 

3n dO{n) 
e" dn " 
and 

r dR{r) T r^ dr{T) 



(10) 



R dr r £2 + r-2 dT 

Equation (10) leads to the solution as; 

e(n) = Ci{n)i 



(11) 



(12) 



3 DESCRIPTION OF CORRELATION 

FUNCTION AND CORRELATION ENERGY 
FOR EXTENDED MASS STRUCTURES 

The solution defined by equation (18) can have variety of 
forms, depending upon the parameters Z and Zn, but we 
are interested in such a solution that is physically valid and 
satisfies earlier results also. A set of boundary conditions 
assumed here are as; 

• The gravitational clustering of galaxies in a homoge- 
neous universe requires ^2 to have a positive value, which 
obviously depends upon the limiting values of n, T, r and e. 

• When nT~^ is very large, the two-particle correlation 
function ^2 will increase and measuring correlation energy 
fee will also increase, and vice versa. 

So depending upon the boundary conditions of ^2 as per 
the requirement, we choose the values of Z and Zn- The 
substitution of equation (18) in equation (3) gives; 



2^Gm^ /3iVy 
3r \4n J 



2 / 2 I 2 \ - 
r (r + e ) 



'dr (19) 



It can be seen, that for different values of Zn, the integral 
will have different values. We fix a set of values for Zn 
starting from unity (1) and check the correlation energy for 
other values also like 0, -1, -2 etc. The following cases are 
discussed here for illustration. 



For Zn = 1, equation (19) leads to; 
27rGm^n ^ /3iV\ # R^ 



3 -V4W 3 (2°) 

Equation (20) relates the variation of (the measuring 
correlation parameter for extended mass structures) with 
the size of cluster. For a given cluster cell with more 
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dimensions of R, we can study the clustering rate without 
involving the thermodynamic quantities. However, there is 
need to show the temperature dependence of be also, as we 
have assumed the system in quasi-equilibrium state. This is 
achieved by testing for other values of Zn- It is interesting 
to note here that the validity of Zn is based on Zn ^ 1- 



For Zn = 0, equation (19) becomes as; 
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For Za 



-1, we have; 
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Finally for Zn ~ —2, we have; 
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The equations (22), (24), and (26) are in good agree- 
ment with a set of boundary conditions. The correlation 
is maximum, if galaxies are treated as point mass objects 
(e = 0) and goes on decreasing as softening parameter (e) 
goes on increasing, thus clustering decreases. It has been 
found that the dispersion in the radial velocity could be up 
to 1000 Kms-^{e.g coma cluster) and is sufficient to throw 
galaxies in the surrounding voids. The cause for decrease 
of correlation for galaxies with extended structures, may be 
due to dispersion in radial velocity. The precise nature of 
such dependence between be and the increase in radial veloc- 
ity can be treated as one of the important problem to study 
in contact to the correlation function studies for galaxy clus- 
ters. Equation (26) is important in the sense that it clari- 
fies the earlier result as well as justifies the signific a nce o f 
be a nd b. In the earlie r work of ISaslaw fc HamiltonI (|l984l l 
and llabal et al.l l|2006l ). it has been understood that b has a 
specific dependence on the combination fiT^'^ . Here in our 
study, it is also true that be has also the specific dependence 
on nT~^, which means that the clustering takes place at 
moderate level. Z may have always positive values because 
of the reason that for large values of N, be increases. 



ena of galaxy clusters. In this paper, an attempt has been 
made to extend the theory of correlation function from point 
mass to extended mass structures by developing a new set 
of equations on the basis of assigned values of the constant 
Zn {Zn=1, 0, -1, -2). The measuring correlation parameter 
be measures the effectiveness of the clustering rate on the 
basis of the two point correlation function ^2 and the re- 
su lts obtained are in good agreement with the earlier work 
of llgbal et al.l l|2006l ). One of the interesting feature in this 
paper is to note how the theoretical values of Zn coincides 
with the observed values of 7 (1.6-1.8) on the same pattern. 
The other interesting result {be having specific dependence 
on the combination nT~^) for the extended structures con- 
firms the earlier results for point mass approximation where 
the functional form of 6 = b{nT~'^). The solution of equation 
(7) given by equation (18) can be used now as an alternative 
tool for studying the phenomena of galaxy clusters on the 
basis of assigned values of Zn- 
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4 DISCUSSION 

The two-point correlation function ^2 for point mass galax- 
ies is the lowest possible tool for understanding the phenom- 
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